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ABSTRACT 
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speed  data  using  five  data  sets  from  areas  where  the  DMSP  satellite  had  crossed  within 
1.1  hours  of  the  GEOSAT  satellite.  Regions  of  ocean  where  the  atmosphere  was  free 
of  cloud  water  and  precipitation  have  a  high  correlation  of  approximately  0.995,  while 
areas  with  precipitation  and  cloud  water  have  very  poor  correlation.  Areas  of  high 
water  vapor  content  do  not  appear  to  adversely  affect  the  high  correlation  between  the 
SSM/I  and  ALTIMETER  measured  wind  speeds. 
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I.  INTRODUCTION 

Understanding  the  environment  has  been  a  problem  mankind  has  dealt  with  for 
ages.  Our  knowledge  of  the  earth  was  initially  limited  to  that  which  we  could  see 
with  our  own  eyes.  Gradually  our  ability  to  sense  our  environment  grew  as  new 
instruments  were  developed.  When  space  satellites  came  onto  the  scene,  a  rich 
resource  of  information  concerning  the  world  around  us  suddenly  opened.  From  visual 
and  infrared  images  of  the  earth,  the  truly  dynamic  nature  of  the  planet  we  live  on 
could  finally  be  appreciated. 

The  years  following  the  first  attempts  at  retrieving  environmental  data  from  space 
have  shown  advances  that  parallel  those  made  in  computer  technology.  Devices  have 
been  implemented  which  have  made  use  of  virtually  all  of  the  electromagnetic  (EM) 
spectrum.  Each  remote  sensing  device  was  designed  for  the  purpose  of  detecting  some 
form  of  radiation  either  emitted  or  reflected  by  the  earth.  To  observe  particular 
features  of  the  atmosphere  or  ocean,  methods  were  developed  which  would  utilize 
particular  frequencies  to  penetrate  through  obstacles  in  the  path  of  this  radiation  and 
arrive  at  its  target  parameter. 

Wind  speed  over  the  ocean  has  proven  to  be  a  particularly  difficult  parameter 
to  measure,  but  most  worthwhile  to  pursue.  Vast  areas  of  ocean  have  been  "data 
deserts"  for  years.  Remote  sensing  of  these  regions  using  microwave  frequencies  can 
supply  data  to  improve  "data  starved"  numerical  input  for  environmental  models. 
This,  in  turn,  should  offer  a  better  understanding  of  synoptic  weather  systems, 
improving  on  the  visible  and  infrared  technologies  currently  used  today.    Disseminating 
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information  concerning  the  global  surface  wind  field  to  all  vessels,  would  enhance 
preparations  for  oceanic  transit.  Ships  already  underway  would  be  better  prepared  for 
heavy  winds  and  seas. 

In  the  case  of  wind  speed  over  the  ocean,  two  separate  techniques  have  been 
developed.  The  first  is  the  active  radar  device  known  as  the  GEOSAT  altimeter  which 
utilizes  pulses  of  EM  energy  at  13.7  GHz  to  analyze  the  surface  of  the  ocean.  Though 
its  prime  objective  was  to  map  the  earth's  geoid,  its  signal  can  also  be  used  to  measure 
surface  wave  height  and  wind  speed  over  the  ocean.  The  second  device  known  as  the 
SPECIAL  SENSOR  MICROWAVE  IMAGER  (SSM/I)  utilizes  a  totally  different 
approach  to  the  problem  of  wind  speed  measurement.  Instead  of  transmitting  its  own 
radiation  to  the  earth  and  then  receiving  a  returned  signal  for  analysis,  it  passively 
measures  upwelling  radiation  from  the  earth's  atmosphere  and  oceans. 

How  these  two  methodologies  compare,  is  the  subject  of  this  thesis.  The  data 
utilized  in  this  comparison  represents  several  areas  of  the  world  where  both  satellites 
were  co-located  within  one  hour  of  each  other's  orbit.  In  the  pages  to  follow,  details 
of  how  each  system  operates  and  the  physics  behind  the  analysis  of  the  data,  will  be 
explored.  Correlation  studies  made  for  each  region  will  disclose  similarities  and 
differences  in  the  wind  speed  measurements.  A  particularly  interesting  case  will  be 
presented  where  several  environmental  parameters  have  interfered  with  the  sensors 
ability  to  accurately  access  the  surface  wind  speed. 

Chapter  two  details  the  basic  physics  applicable  to  the  SSM/I  and  GEOSAT 
remote  sensors,  as  well  as  the  operation  of  both  systems.  Chapter  three  will  discuss 
ground  truth  verification  studies  made  for  each  system.    Here,  individual  strengths  and 


weaknesses  are  taken  into  account  prior  to  the  comparison  study.  The  experiment  is 
presented  in  chapter  four,  followed  by  conclusions  and  recommendations  in  chapter 
five. 


II.  FUNDAMENTALS  OF  MICROWAVE  REMOTE  SENSING 

The  basic  physics  behind  the  remote  sensing  techniques  examined  in  this  thesis, 
will  be  discussed  in  this  chapter.  The  first  section  will  cover  the  propagation 
characteristics  of  EM  radiation.  This  will  be  followed  by  an  overview  of  the  index 
of  refraction,  blackbody  radiation,  emissivity,  and  the  radiative-transfer  equation.  This 
should  provide  an  ample  background  for  understanding  the  operation  of  both  the  SSM/I 
and  the  GEOSAT  altimeter. 

A.   PROPERTIES  OF  EM  WAVE  PROPAGATION  AND  ATTENUATION 

Microwave  remote  sensing  depends  on  the  modification  of  EM  radiation  by  the 
environment.  Changes  take  place  in  the  phase  and  intensity  of  an  EM  wave  as  it 
passes  through  the  atmosphere  and  interacts  with  the  ocean  surface.  These  alterations 
of  the  wave's  character  are  interpreted  by  the  remote  sensing  device  as  variations  in 
the  state  of  the  environment.  We  will  first  describe  the  EM  wave  equation  and  then 
develop  the  relationships  that  govern  the  changes  that  take  place  in  the  wave  as  it 
passes  through  matter. 

EM  radiation  consists  of  alternating  electric  and  magnetic  fields  [Ref.  1],  which 
are  classified  by  their  frequency  of  fluctuation.  A  single  EM  wave  vector,  E(x.y,z), 
can  be  broken  down  into  its  individual  components,  EK,  Ey,  and  Ez  for  a  simplified 
analysis  (Fig.l).  The  propagation  characteristics  of  a  single  component  of  an  EM  wave 


Direction  of 
propagation 


Figure  1  Electromagnetic  wave  vector  and  components. 


can  be  described  by  the  following  equations: 


E   =  E   eH(fa-",J 
k2  =  co2u£  -  icou 
co  =  2rtf 

e  =  £0er  =  £„(£'  -  ie") 


where: 


Ey  s  component  of  the  electric  field  intensity 

Ec  =  initial  electric  field  magnitude 

k  =  complex  wave  number 

co  =  frequency  in  radians  per  second 

f  =  frequency  in  hertz 


(1) 

(2) 
(3) 
(4) 
(5) 


(i0  =  permeability  of  free  space  (4k  x  10'7Hz/m) 

|Xr  =  relative  permeability  of  a  material  (=1) 

e  =  permittivity  of  a  medium 

e0  =  permittivity  of  free  space  (8.85  x  10'12  Farads/m) 

e',£"  =  real  and  imarginary  components  of  the  dielectric  constant 

For  most  cases  in  remote  sensing  applications,  the  value  of  fir    =  1,  so  that 
equation  (2)  can  be  written  as: 

k2  =  co2u0e0[e'  -  ia/(eca»]  (6) 

where: 

er  =  [e'  -  ia/(eo(o)]  the  complex  dielectric  constant 

G  =  conductivity  of  the  medium 

In  free  space,  i.e.,  a  vacuum,  e,  =  1,  a  =  0,  and  the  EM  wave  velocity  is  given 
as: 

c02  =  l/(u0e0)  (7) 

For  EM  waves  traveling  through  matter,  such  as  the  ocean  or  atmosphere,  the  velocity 
is  determined  through  the  equation: 

c  =  co/k  =  cjJz,  (8) 

From  (7)  and  (8),  we  can  define  the  index  of  refraction  n,  which  is  related  to  er  by 
the  formula: 

n  =  cjc  =  V£  (9) 

Equation  (9)  is  important  for  understanding  how  the  EM  wave  is  effected  as  it  passes 
through  the  atmosphere  and  interacts  with  the  ocean  below.  Since  n  is  related  to  the 
complex   dielectric   constant  er,   radiation   passing  through  matter  will   experience   a 


modification  in  speed  and  intensity.  The  wave  will  be  slowed  due  to  the  resistance 
encountered  by  an  induced  electric  field  produced  in  this  medium  and  attenuated 
through  the  imaginary  component  of  n.  To  understand  how  the  imaginary  component 
can  attenuate  the  wave,  n  must  be  examined  in  more  detail.  The  index  of  refraction 
written  in  its  component  form: 

n  =  a  +  ip  (10) 

where  a  and  p  represent  the  real  and  imaginary  parts  of  n  and  are  written  as: 

a2  =  (|er|  +  e')/2  (11) 

P2  =  (|er|  -e')/2  (12) 

IeJ  =  [e,2  +  c7(e0©)T  (13) 

Returning  to  equation  (8)  for  the  wave  number  k  and  substituting  the  component  form 
of  n  into  that  expression,  we  can  rewrite  k  as: 

k  =  uxx/c    +  icop/c  (14) 

By  replacing  the  expression  for  k  in  equation  (14)  into  the  original  expression  for  the 
propagating  wave,  equation  (1),  we  have: 

£      _    £      e[-i((oap/c     -tot)]    e[«Px/c   ]  (15) 

_    C       gl**]     e[-l(con*/c     -  rot)J 

where  ({)  =  -u)/c  is  known  as  the  attenuation  coefficient.  With  a  negative  value  in  the 
exponential,  this  component  acts  to  reduce  the  value  of  the  electric  field  as  the  EM 
wave  passes  through  matter.  The  distance  the  radiation  travels  through  some  medium, 
in  which  the  intensity  of  the  wave  is  reduced  by  a  factor  of  exp(-l),  is  known  as  the 
skin  depth  d  and  is  given  by  : 

d  =  c/((0p)      or      d  =  c/(2jcfp)  (16) 


For  a  frequency  of  5  GHz,  the  value  of  d  is  0.5  centimeters.  The  SSM/I  operates 
between  19  and  37  GHz  and  the  GEOSAT  altimeter  operates  at  13.5  GHz,  and  so  the 
skin  depths  for  the  instruments  in  this  study  are  on  the  order  of  millimeters. 

B.       POLARIZATION 

In  addition  to  the  variation  of  the  intensity  and  frequency  of  the  wave,  microwave 
radiation  is  also  characterized  by  its  polarization.  Polarization  is  defined  as  the 
orientation  of  the  electric  field  vector  relative  to  the  surface  it  is  reflected  or 
emitted  [Ref.  2].  Polarized  radiation  is  classified  as  either  horizontal,  in  which  case 
the  field  is  parallel  to  the  surface,  or  vertical,  whereby  the  field  is  perpendicular  with 
respect  to  the  surface.  A  change  in  surface  conditions  and  viewing  angle,  will  affect 
the  vertically  and  horizontally  polarized  EM  waves  in  different  ways.  The  amount  of 
radiation  emitted  from  a  surface  (emissivity)  depends  on  how  the  radiation  is  polarized 
(Fig. 2).  The  SSM/I's  sensors  can  detect  both  horizontally  and  vertically  polarized 
radiation.  Certain  combinations  of  polarized  frequencies  are  used  to  allow  particular 
environmental  parameters  to  be  measured  separately  with  a  minimal  of  interference 
from  the  other  parameters. 

Assuming  a  smooth     (specular)  surface,  the  reflection  of  the  spectral  radiance, 
R,  from  such  a  surface  is  written  in  equation  form  as: 

R(A.;er,4>r-7c)  =  Lra;er,<t>,)/LA;e,<K-K)  U7) 

where: 

X  =  wavelength  of  the  radiation 
9,<}>  =  spherical  coordinate  angles 
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Figure  2    Polarized    microwave    radiation    at    3    cm    wavelength    as    a 
function  of  emissivity  and  incidence  angle. 


Lr,L,  =  reflected  and  incident  radiance  respectively 


The  spectra]  reflectance  of  a  plane  dielectric  surface  is  determined  through  the  Fresnel 
reflection  coefficients  as  follows  [Ref.  3]: 

RH  =  [(p  -  cos6)2  +  q2]/[(p  +  cos0)2  +  q21  (18) 

Rv  =  [(e'cosG  -  p)2  +  (e"cosG  +  q)2]/[(e'cos0+  p)2  +  (e"cos9  +  q)2]         (19) 

where  Rv  and  RH,  are  the  vertically  and  horizontally  polarized  components  of  the  EM 

wave,  and  e  =  e'  -  ie";  the  complex  dielectric  constant  of  the  surface.  The  quantities 

p  and  q  are: 

p  =  (1/V2){[(£'  -  sin29)2  +E"T  +  [£'  -  sin26]p  (20) 


p  =  (1/V2){[(e'  -  sin20)2  +e"T  +  [e'  -  sin20]P  (20) 

q  =  (1/V2){[(e'  -  sin20)2  +e"2f  -  [e'  -  sin2©]}*  (21) 

where  6  is  the  incidence  angle. 

In  the  case  when  0  =  0,  such  as  with  the  viewing  angle  of  the  GEOSAT 
altimeter,  RH  and  Rv  are  reduced  to  the  form: 

RH  =  Rv  =  [(p  -  I)2  +  qV[(p  +  D2  +  q2]  (22) 

=  [(a  -  l)2  +  (32/[(a  +  l)2  +  p2] 
where  the  complex  index  of  refraction  is  given  by  equation  (10).  Since  n2  =  e'  -  ie", 
the  Fresnel  reflection  coefficient  at  nonnal  incidence      can  take  the  form: 

Rv  =    |(n  -  l)/(n  +  1)|2  (23) 

The  spectral  emissivity  of  a  surface  is  related  to  the  spectral  reflectivity  through 
Kirchoff's  law  (20).  We  can  relate  the  spectral  reflectance  R  to  the  emissivity  in 
both  horizontally  and  vertically  polarized  forms: 

ev(6)  =  1  -  Rv(9)  (24) 

e„(e)  =  1  -  RH(6)  (25) 

C.       BRIGHTNESS  TEMPERATURE 

The  purpose  of  this  section  is  to  develop  the  relationship  between  the  radiation 
emitted  by  the  ocean's  surface  and  the  upwelling  temperature  detected  by  the  SSM/I 
known  as  the  brightness  temperature,  Tb.  All  objects,  except  those  existing  at  absolute 
zero,  emit  radiation.  An  object  is  classified  as  a  perfect  emitter,  known  as  a  black 
body,  if  it  radiates  energy  at  the  maximum  rate  predicted  by  Planck's  radiation  law. 
As  the  temperature  of  an  object  increases,  the  total  amount  of  energy  emitted  also 
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increases.  Most  objects,  however,  emit  at  a  rate  less  than  that  of  a  blackbody.  The 
term  "grey  body",  is  used  to  describe  an  object  which  emits  radiation  at  a  reduced  rate. 
Variations  in  the  temperature  of  a  grey  body  relative  to  the  assumed  blackbody 
temperature  of  the  earth  can  be  identified  and  attributed  to  changes  in  the  environment, 
such  as  wind  speed,  water  vapor,  and  precipitation. 

Terms  commonly  used  in  discussing  brightness  temperature,  or  radiation 
measurement  in  general,  are  radiant  energy,  radiant  flux,  irradiance,  and  emittance. 
The  term  radiant  energy,  Q,  refers  to  the  amount  of  energy  carried  by  the  radiation 
(Joules).  The  radiant  flux,  O  =  dQ/dt,  is  a  measure  of  rate  of  energy  transport 
(Joules/sec).  This  flux  defines  the  power  available  to  produce  a  response  in  a 
measuring  device.  The  irradiance,  E  =  dO/dm2,  is  used  in  reference  to  the  radiance 
incident  on  a  surface.  The  radiance  emitted  from  a  surface  is  known  as  the  emittance, 
M  =  dO/dt.  The  character  of  the  material  from  which  radiation  is  either  emitted  or 
reflected  is  known  as  its  emissivity.  Each  surface  has  a  particular  value,  for  example, 
the  emissivity  of  salt  water  is  0.3,  while  that  of  dry  soil  is  0.92.  It  is  this  difference 
in  emissivity  that  allows  the  SSM/I  to  distinguish  between  land  and  sea  surfaces. 

The  GEOSAT  relies  on  the  degree  of  reflectivity  to  detect  changes  in  the  ocean's 
surface.  The  GEOSAT  and  SSM/I  observations  are  linked  by  Kirchoff  s  law,  which 
relates  emissivity  and  reflectivity  in  the  following  manner: 

e  =  1  -  r  (26) 

where  r  is  the  reflectivity  of  the  surface  and  e  is  the  emissivity. 

The  ratio  of  spectral  brightness  emitted  in  a  direction  (0,<|))  by  a  real  body 
relative  to  that  which  is  radiated  by  a  blackbody  is  the  spectral  directional  emissivity 
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e  of  that  body,  where  E  is  a  function  of  wavelength  X  and  spherical  angles  (j)  and  0. 
The  microwave  emission  from  the  ocean  surface  is  given  by  eBv,  where  Bv  follows 
the  Rayleigh- Jeans  approximation  of  Planck's  Law.  Planck's  Law  states  that  the  spectral 
brightness  per  unit  of  frequency  is  expressed  as: 

Bv  =  (2hv2/c2)  (I/re'™"™  -  1])  (27) 

where: 

k  =  1.38  x  10"23  J/K,  Boltzman's  constant 

h  =  6.63  x  10"34  J  s,  Planck's  constant 

T  =  temperature  in  degrees  Kelvin 

v  =  frequency  of  radiation  in  hertz 

c  =  3  x  108  m/s,  speed  of  light  in  a  vacuum 
The  Rayleigh-Jeans  approximation  is  made  at  microwave  frequencies,  where  hv/(kT) 
«  1,  so  that: 

e.hv/(kT„     =     J     +    hv/(kT) 

Substituting  the  above  relationship  into  equation  [18]  gives 

Bv  =  (2v2k/c2)T  (28) 

Equation  [19]  states  that  the  spectral  brightness  is  proportional  to  the  temperature  of 
the  body  emitting  radiation.  The  radiance  received  by  a  microwave  sensor  is  known 
as  the  brightness  temperature.  The  brightness  temperature  Tb,  is  given  by  eT  and  will 
be  discussed  in  a  later  chapter. 

The  following  section  will  cover  the  topics  of  absorption  and  the  interaction  of 
radiation  with  surfaces.  These  factors  contribute  to  the  variation  in  brightness 
temperatures  as  measured  from  satellites. 
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I).       ABSORPTION  AND  SCATTERING  OF  RADIATION 

Radiation  emitted  by  the  ocean  surface  must  pass  through  the  atmosphere  once 
before  being  detected  by  the  SSM/I,  and  twice  for  the  GEOSAT  altimeter. 
Consequently,  interactions  which  take  place  between  the  atmosphere  and  the  radiation 
must  be  understood  in  order  to  interpret  a  retrieved  signal.  In  the  atmosphere,  energy 
is  both  scattered  and  absorbed.  For  microwave  frequencies,  the  primary  absorber  is 
water  vapor,  while  scattering  is  predominantly  due  to  rain.  Any  instrument  used  to 
measure  wind  speed  over  the  ocean  must  utilize  frequencies  where  the  effects  of 
atmospheric  absorption  and  scattering  are  small.  The  GEOSAT  altimeter  frequency  was 
chosen  such  that  atmospheric  effects  would  be  small.  In  addition  to  observing  the 
ocean  surface,  the  SSM/I  measures  atmospheric  parameters  such  as  water  vapor  and 
precipitation.  This  allows  the  SSM/I  to  adjust  for  atmospheric  parameters  which  would 
interfere  with  the  measurement  of  wind  speed.  Questionable  wind  speeds,  based  on  the 
presence  of  these  interfering  parameters,  are  identified  at  ground  stations  when  the  data 
is  analyzed. 

E.       THE  RELATION  OF  SURFACE  ROUGHNESS  TO  EMITTED  RADIATION 

The  roughness  of  the  ocean  causes  the  reflected  and  emitted  EM  waves  to  be 
sent  out  in  numerous  directions.  Remote  sensors  of  microwave  radiation  must  be  able 
to  account  for  this  variation  in  directionality,  hence  the  viewing  angle  becomes  an 
important  factor.    The  received  radiation  is  also  a  function  of  the  dielectric  constant 


13 


of  the  surface.  These  properties  can  be  related  if  we,  (1)  assume  that  the  surface  is 
made  up  of  many  small  facets,  (2)  examine  only  a  single  facet,  or  plane  dielectric 
surface  and  apply  the  Fresnel  Reflection  Coefficients  to  that  surface  to  obtain  the 
spectral  reflectance  of  that  surface.  Finally,  (3)  integrate  over  all  the  facets  in  the  area 
observed  by  the  sensors  to  obtain  the  spectral  reflectance  for  that  surface. 

The  solution  to  the  Fresnel  equations  requires  that  either  the  complex  dielectric 
constant,  er,  or  the  index  of  refraction  of  the  surface  as  a  function  of  frequency  be 
known.  In  the  case  of  the  ocean,  the  dielectric  constant  is  a  well  known  parameter 
through  the  work  of  Lane  and  Saxon  [Ref.  4].  The  dielectric  constant  is  represented 
by  the  Debye  equation: 

er  =  e_  +  [(£.£j/(l  +  (ioox)1^  -  ia/ooej  (29) 

where: 

(0  =  27cf,  radiation  frequency 

£„,  =  dielectric  constant  at  infinite  frequency 

£,  =  the  static  dielectric  constant 

x  =  the  relaxation  time 

a  =  the  ionic  conductivity 

a  =  an  empirical  constant 
The  parameters  £„  T,  and  a  are  all  functions  of  the  material  which  is  reflecting  or 
emitting  the  radiation  to  be  received  by  microwave  sensors.     The  conductivity  o  of 
the  sea  water  increases  the  complex  portion  of  the  dielectric  constant,  however,  it  is 
important  to  note  that  the  influence  of  sea  salts  on  the  emissiviry  of  the  ocean  is  small 
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for  frequencies  greater  than  5  GHZ.     Since  both  the  SSM/I  and  GEOSAT  sensors 
operate  above  5  GHZ,  the  salinity  of  the  sea  water  has  no  effect. 

The  above  description  for  spectral  reflectance  and  emittance  from  a  single  facet 
can  now  be  applied  to  the  entire  rough  surface.  This  is  accomplished  by  integrating 
over  all  the  facets  and  taking  a  statistical  average. 

F.       THE  RADIATIVE  TRANSFER  EQUATION 

The  radiative  transfer  equation  contains  terms  which  relate  radiant  energy  from 
the  surface  and  atmosphere  to  the  radiant  energy  observed  by  a  satellite.  Before 
energy  can  reach  a  satellite  sensor  from  the  ocean  surface,  the  upwelling  radiation  is 
attenuated  and  joined  by  additional  emission  from  the  atmosphere,  so  that  the  observed 
temperature  detected  by  the  SSM/I  is  given  by: 

T  =  TA(1  -  t)  +  t[eT.  +(1  -  e)(l  -  t)TJ  (30a) 

or 

T  =  etfT.  -  TA(l-t)]  +  TA(1  -  t2)  (30b) 

where,  TA  is  the  effective  temperature  of  the  atmosphere,  t  =  e'"TI,  T  is  the  optical 
depth,  T,  is  the  sea  surface  temperature,  and  e  is  the  emissivity.  The  brightness 
temperature  of  the  radiation  incident  on  the  ocean  surface  by  the  atmosphere  is 
TA(  It).  The  amount  that  is  reflected  by  the  sea  surface  is  written  as  (1  -  e)TA(l  - 
t),  while  the  amount  emitted  is  eT,.  It  follows  that  the  total  upwelling  brightness 
temperature  from  the  sea  surface  is  eT,  +  (1  -  e)(l  -  t)TA.  The  signal  from  the 
surface  is  ET,.     The  other  terms  are  considered  noise  with  regard  to  surface  wind 
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speed  analysis,  but  are  useful  when  detecting  water  vapor,  liquid  water,  and 
precipitation.  If  the  optical  depth  x  =  0,  then  the  atmosphere  is  transparent  and  allows 
just  the  brightness  temperature  of  the  ocean  surface  to  be  seen  by  the  sensor.  Equation 
(30b)  would  then  be  reduced  to  simply  T  =  eT,.  This  the  ideal  case  for  measuring  wind 
speed  over  the  ocean.  If,  on  the  other  hand,  x  »  0,  then  T  =  TA,  i.e.,  the  SSM/I 
would  be  looking  at  the  environmental  changes  in  the  atmosphere  only.  The 
frequencies  and  the  method  used  by  the  SSMI  and  GEOSAT  to  minimize  atmospheric 
interference  will  be  discussed  in  sections  devoted  to  instrument  operation. 

G.  THE  TRANSFER  OF  MOMENTUM  FROM  THE  ATMOSPHERIC  SURFACE 
LAYER  TO  THE  SEA  SURFACE 

Waves  and  foam  generated  by  the  wind  contribute  to  the  roughness  of  the  sea 
surface  and  effect  both  microwave  emission  as  measured  by  the  SSM/I  and  microwave 
scattering  as  measured  by  the  GEOSAT  altimeter.  The  roughness  has  the  effect  of 
altering  the  incidence  angle  of  radiation  being  emitted  or  reflected  and  changing  the 
emissivity  of  the  sea  surface  itself.  The  roughness  of  the  sea  surface  is  an  indication 
of  the  transfer  of  momentum  into  the  sea  from  the  turbulent  atmospheric  surface  layer. 
An  extrapolation  process  is  performed  to  relate  the  wind  speed  in  the  atmospheric 
surface  layer  to  the  sea  surface  roughness.  This  section  will  briefly  describe  the 
transfer  of  momentum  from  the  atmospheric  surface  layer  to  the  sea  surface. 

Businger  [Ref.  5]  describes  the  formulation  of  the  turbulent  transfer  of  momentum 
from  the  atmospheric  surface  layer  to  the  sea  surface  as  a  semi-empirical  relationship. 
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The  appropriate  equation  can  be  developed  by  utilizing  the  turbulent  transfer  coefficient 
K  in  the  equation  for  momentum  transfer,  where: 

Fm  =  -pKaTJ/az  (31) 

Here,  p  is  the  density  and  U  is  the  mean  horizontal  wind  speed.  Businger  defines 
the  atmospheric  surface  layer  as  the  lower  20  meters  of  the  atmosphere.  Here  the 
interaction  with  the  surface  is  strong  and  the  layer  adjusts  to  the  surface  conditions 
rapidly.  We  can  therefore  assume  a  quasi-steady-state.  The  momentum  transferred 
through  to  the  ocean  surface  exerts  a  drag  force  per  unit  area  on  that  surface  known 
as  the  surface  stress  x.    By  Newton's  third  law,  the  momentum  equation  is  given  as: 


Fm  =  puV  =  -T  (32) 
The  friction  velocity  u.  is  commonly  used  as  a  scaling  velocity,  where: 

u.  =  (x/p)m  (33) 
Equation  (31)  can  now  be  written  as: 

KjU/dz  =  u.2  (34) 

where  u.  is  independent  of  height.    To  generate  a  wind  profile  from  the  surface  to  a 
height  of  20  meters,  an  assumption  must  be  made  about  K^.    By  dimensional  analysis, 

K,,,  has  the  dimensions  of  velocity  times  length.    We  use  u.  to  scale  the  velocity  and 
z  to  scale  the  length.    Therefore,  we  assume  that: 

K.  =  ku.z  (35) 
where,  k  is  a  constant  of  proportionality,  the  von  Karmann  constant.  Combining 
equations  (34)  and  (35),  we  have: 

dU/9z  =  u./(kz)  (36) 
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To  avoid  an  infinite  shear  at  the  surface,  when  z  =  0,  a  finite  surface  roughness  z0  is 
introduced  into  the  equation.  This  now  becomes: 

dU/dz  =  u./(k[z  +  zj)  (37) 

Upon  integration  we  have  the  form: 

U/u.  =  1/k  ln[(z  +  z0)/zj  (38) 

which  is  the  logarithmic  profile  for  the  wind  from  the  sea  surface  through  the 
atmospheric  surface  in  a  neutrally  stable  atmosphere.  This  is  the  common  profile  used 
to  extrapolate  surface  wind  speed  to  a  height  of  10  or  20  meters,  which  are  the  typical 
heights  at  which  sea  buoys  and  ships  measure  wind  speed  at  sea.  The  neutral  case  is 
used  primarily  for  its  simplicity  and  also  for  its  adequate  approximation  of  the  surface 
layer  wind  field.  Figure  (3)  [Ref.  6]  represents  a  semi-log  plot  of  wind  speed  from 
the  surface  through  the  atmospheric  boundary  layer  for  neutral,  stable  and  unstable 
conditions.  It  would  be  impractical  to  measure  the  stability  of  the  atmosphere  for  all 
wind  speed  measurements.  Fortunately,  the  neutral  profile  of  wind  speed  approximates 
the  atmosphere  in  most  situations. 

Fig.  3  is  derived  from  setting  equation  (37)  equal  to  a  wind  shear  function,  Om, 
such  that:  <X>M  =  (k[z  +  zj/u.)  dU/dz 

where: 

Om  =  1  +  [  4.7  z/L  ],  for  z/L  >  0  (stable)  (39a) 

Om  =  1  ,  for  z/L  =  0  (neutral)  (39b) 

<DM  =  [1  +  [  15  z/L  ]]-,/4,        for  z/L  <  0  (unstable)  (39c) 

and  L  is  the  Obukhov  length. 
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Figure  3    Typical    wind    speed    profiles    versus    static 
stability  in  the  surface  layer. 
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III.  MEASURING  WIND  SPEED  OVER  THE  OCEAN  WITH  THE  SSM/I 

The  SSM/I  is  a  microwave  radiometer  which  measures  upwelling  radiation  from 
the  earth's  surface  at  five  frequencies  and  two  polarizations.  Using  frequencies 
associated  with  atmospheric  windows,  the  SSM/I  measures  upwelling  thermal 
microwave  radiation  from  the  ocean  surface.  This  radiation  is  called  the  brightness 
temperature.  From  equation  (30b),  where  Tb  =  eT,  for  an  atmosphere  with  100% 
transmittance,  the  brightness  temperature  is  a  function  of  the  emissivity,  e,  and 
temperature  of  the  sea  surface,  T,.  The  sea  surface  temperature  is  relatively  constant 
compared  to  the  change  in  emissivity.  An  increase  in  Tb  is  thus  caused  by  an  increase 
in  the  sea  surface  emissivity.  A  change  in  surface  roughness  and  thus  wind  speed,  will 
change  the  emissivity  of  the  surface.  Thus,  a  particular  wind  speed  at  the  sea  surface 
can  be  related  to  a  certain  brightness  temperature. 

The  SSM/I  utilizes  the  19.35H,  22.235V,  37H,  and  37V  GHz  channels  for 
measuring  wind  speed,  where  H  and  V  represent  horizontal  and  vertical  polarizations 
respectively.  The  19.35H,  37H,  and  37V  GHz  channels  are  "windows"  in  which  there 
is  nearly  100  %  transmission  of  upwelling  radiation  through  the  atmosphere.  The 
22.235  GHz  channel  is  in  the  center  of  a  water  vapor  absorption  line  and  is  used  to 
measure  water  vapor  concentration.  Detection  of  water  vapor  is  vital  to  the 
determination  of  the  validity  of  the  wind  speed  detected,  since  the  19.35H,  37V,  and 
37H  channels  are  subject  to  some  small  amount  of  atmospheric  interference  which 
require  that  an  atmospheric  correction  be  applied. 


20 


The  algorithm  used  by  the  SSM/I  for  measuring  wind  speed  assumes  the 
loughness  of  the  ocean  surface  can  be  treated  as  a  conglomeration  of  many  tiny 
emitting  planes  whose  angle  of  inclination  is  altered  by  the  force  of  the  wind.  The 
variation  in  the  sea  surface  slope  is  used  in  calculating  the  emissivity  of  the  surface 
by  assuming  a  Gaussian  slope  distribution  of  the  angles  of  the  slope.  A  change  in 
the  distribution  of  these  angles  corresponds  to  a  change  in  the  root  mean  square  slope 
of  the  sea  surface  and  distribution  of  the  emissivities  associated  with  the  variance  of 
sea  surface  slopes.  The  resulting  change  in  emissivity  causes  a  change  in  brightness 
temperature  through  the  Fresnel  equations  (Eqns.22  and  23).  The  equation  used  to 
describe  the  dependence  of  derived  wind  speed  on  the  surface  slope  variance,  & 
[Ref.  7]  ,  is  as  follows: 

<j2(f)  =  (0.3  +  0.02f)  (0.003  +  0.48U2O)    for  f  <  35  GHz  (40a) 

&(f)  =  0.003  +  0.48U2o    for  f  >  35  GHz  (40b) 

Here,  f  is  defined  as  the  frequency  in  GHz  and  V^  as  the  wind  speed  measured  in 
meters/second  at  a  20  meter  height. 

Foam  generated  by  breaking  waves  occurs  at  wind  speeds  greater  than  7  m/s. 
The  affect  of  foam  on  the  microwave  emission  from  the  sea  surface  is  not  well 
understood  and  there  are  varying  theories  in  literature.  Swift  TRef.  81  cites  an 
experiment  over  the  Salton  sea  using  a  19  GHz  radiometer  to  measure  wind  speed 
over  a  foam  covered  sea.  Swift  states  that  the  results  were  tens  of  degrees  higher  than 
theory  and  could  be  attributed  to  the  presence  of  foam.    Wilheit  [Ref.  9]  treated  foam 
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as  a  relatively  non-reflective  surface  and  eliminated  it  from  the  radiative  transfer 
equation.  The  reflectivity  that  does  occur  comes  from  the  dielectric  interface  around 
the  foam  bubbles.  Hollinger  et  al  [Ref.  10]  treat  foam  as  "a  substance  which  partially 
obscures  the  ocean  surface  in  a  manner  independent  of  polarization  but  dependent  upon 
frequency".  He  defines  the  foam  fraction,  Ff,  as 

Ff  =  0.006  (Ujo  -  7)  (1  -  e(n-5  )      for  UM  >  7  m/sec  (41a) 

Ff  =  0    for  U^  <  7  m/sec  (41b) 

This  equation  is  developed  from  the  observation  that  the  sea  surface  brightness 
temperature  increases  linearly  for  wind  speeds  in  excess  of  7  m/sec,  and  that  there  is 
an  absence  of  foam  development  for  speeds  less  than  or  equal  to  7  m/s.  The  foam 
fraction  is  also  used  in  the  empirically  derived  SSM/I  wind  speed  algorithm  to  account 
for  the  change  in  emissivity  that  occurs  with  the  presence  of  foam. 
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IV.  MEASURING  WIND  SPEED  OVER  THE  OCEAN  WITH  THE  GEOSAT  ALTIMETER 


The  primary  mission  of  the  GEOSAT  satellite  was  not  wind  measurement,  but 
measurement  of  the  geoid  and  sea  surface  height.  However,  wind  is  computed  as  a 
by-product  and  with  an  accuracy  of  +/-  2  m/s  [Ref.  10].  To  measure  wind  speed,  the 
GEOSAT  satellite  utilizes  an  active  radar  operating  at  a  frequency  of  13.9  GHz  called 
an  altimeter.  This  instrument  transmits  a  pulse  of  microwave  energy  to  the  ocean 
surface  and  measures  the  return  signal.  The  altimeter  [Ref.  11]  uses  an  automatic  gain 
control  (AGC)  loop  to  normalize  the  amplitude  of  the  return  signal  from  the  sea 
surface.  The  AGC  is  a  measure  of  the  backscatter  cross  section  of  the  sea  surface,  a  , 
and  is  related  to  the  wind  speed.  The  greater  the  wind  speed,  the  lower  the  cross 
section,  thus  the  weaker  the  return  signal  from  the  ocean  surface.  The  backscatter 
cross  section  is  obtained  from  the  following  equations  [Ref.  12]: 

AGC  [dBm]  =  10  In  (P./PJ  (43) 

where: 

Pr  =  the  received  power 

P„  =  the  reference  power  (1  mW) 

The  received  power  Pr,  is  defined  in  terms  of  antenna  gain  G,  transmitted  power  P„ 
the  operating  wavelength  A.,  the  distance  to  the  sea  surface  R,  and  the  scattering  cross 
section  a  .  The  relation  appears  as  follows: 

pr  =  [\7/(4n)]      P,  G  a  /  R4dA  (44) 
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Assuming  that  the  sea  surface  slope  is  both  gaussian  and  isotropic  in  its  probability 
distribution,  we  can  define  a    as  follows: 

a    =  [  |R(0)I2  /  S2  ]  sec40  exp  [  -tan20  /  S2  ]  (45) 

where: 

0  =  angle  of  incidence 

R(0)  =  Fresnel  reflection  coefficient  for  the  air-sea  interface  at  normal  incidence 

The  Cox   and  Munk   [Ref.  13]  relationship  for  average  upwind/crosswind  mean 
square  sea  surface  slope  a2  and  wind  speed  is  given  as: 

a2  =  a  +  6  U  (46) 

where,  U  is  the  wind  speed  in  m/s  and  a,  6  are  empirical  constants  of  proportionality 
which  depend  on  the  probability  density  of  the  sea  surface  slope.  By  making  the 
assumption  of  an  isotropic  Gaussian  surface,  the  separate  effects  of  the  upwind 
component  and  the  cross  wind  component  are  ignored  [Ref.  14].  By  ignoring  wind 
direction  when  deriving  a  relationship  between  wind  speed  and  surface  slope,  error  in 
the  surface  wind  speed  will  be  introduced. 

If  we  let  0  =  0  ,  since  we  are  observing  at  the  nadir  angle,  the  equation  for  a 
becomes: 

a   =  lR(0)l2/(a  +  BU)  (47) 

The  equation  for  the  received  power  is 

Pr  =  (X2/(4tc)3)  P,  (G7h4)  a   A  (48) 

where: 

h  =  satellite  altitude 
A  =  area  illuminated  by  the  altimeter 
This  relationship  is  used  in  the  Brown  algorithm  for  the  actual  calculation  of  wind 
speed  [Ref.  15]  at  the  Fleet  Numerical  Oceanographic  Center.    Brown's  algorithm  is 
a  result  of  a  comparison  of  wind  speed  data  received  from  sea  buoys  to  normalized 
backscattered   cross   sections  derived   from   altimeter  power  measured   on  board   the 
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GEOS-3  satellite  and  applied  to  GEOSAT.  The  algorithm  is  presented  in  two  stages. 
The  first  section  is  as  follows: 

W,  =  exp  [(io-<02l*«™°>  -  B)/A]  (49) 

where: 

A  =  0.080074    B  =  -0.124651    for  a  <  10.12  dB 

A  =  0.039893    B  =  -0.031996    for  10.12  <  a*  <  10.90  dB 

A  =  0.015950    B  =    0.017215    for  a'  >  10.9  dB 

In  the  above  equation,  the  backscattered  cross  section  is  given  as  decibels  relative  to 

1  m2. 

A  second  equation  is  used  to  complete  the  derivation  of  surface  wind  speed.  This 

was  necessary  to  adjust  for  the  skewness  in  the  distribution  of  the  differences  between 

the  buoy  and  altimeter  wind  speed  measurements.  The  second  part  of  the  algorithm  is 

as  follows: 


U10  =  Ia,W, 

n-1 


for  W,  <  16  meters  per  second 
for  W,  >  16  meters  per  second 


(50a) 

(50b) 


where: 


a,  =    2.08779900000 

a2  =  -0.36499280000 

a3  =    0.04062421000 

a4  =  -0.00190495200 

a,  =    0.00003288189 

U,0  is  the  wind  speed  at  a  height  of  10  meters.     This  algorithm  appears  to  measure 
the  wind  to  within  +/-  2  m/s. 
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V.  GROUND  TRUTH  STUDIES 

A.  GEOSAT  Wind  Speed  Ground  Truth  Study 

Shuhy  [Ref.  16]  made  a  ground  truth  study  for  the  wind  speeds  obtained  by 

GEOSAT  using  National  Data  Buoy  Center  buoy  data.  To  avoid  possible  interference 
from  land,  most  of  the  buoys  used  were  located  in  deep  water  away  from  any 
coastlines.  Shuhy  chose  to  restrict  the  comparison  to  satellite  passes  within  80 
kilometers  and  90  minutes  of  a  ground  truth  point.  The  data  was  classified  as  either 
excellent,  good,  or  fair  based  on  the  following  criteria: 


DISTANCE  FROM 

TIME  DIFFERENCE 

CLASSIFICATION 

GROUND  TRUTH 

IN  OBSERVATIONS 

excellent 

25  KM 

30  min 

good 

25  -  60  KM 

60  min 

fair 

60  -  80  KM 

90  min 

No  data  outside  80  kilometers  or  90  minutes  from  crossing  were  used.  The  results  of 
the  comparison  are  seen  in  Fig.  4.  Shuhy  concluded  from  regression  analysis  that  the 
sub-categorization  described  above  yielded  no  significant  differences.  A  total  of  256 
data  points  were  used  with  an  rms  error  of  1.54  m/s.  The  comparison  demonstrated  that 
GEOSAT  wind  speeds  are  higher  than  the  ground  truth  for  speeds  less  than  6  m/s  and 
for  speeds  greater  than  6  m/s,  GEOSAT  winds  are  lower  than  those  obtained  on  the 
surface.  Shuhy  notes,  however,  that  there  is  significant  scatter  in  the  data  and  that 
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■  ▲•  Regression  line 
■  *    I     «    ■ 


■    0  to  25  km,  30  min 

*  25  to  60  km,  60  min 

•  60  to  80  km,  90  min 


3  6  9  12  15 

Ground-truth  wind  (meters  per  second^ 
4  Geosat  vs  ground  truth  winds  at  a  height  of  10  meters. 


spatial  variability  must  be  examined  to  determine  which  comparisons  might  be  invalid 
due  to  high  variability  between  buoy  and  satellite  wind  speeds  along  the  satellite  track. 
When  ground  truth  data  is  compared  to  off-nadir  GEOSAT  data,  there  appears  to 
be  a  correlation  between  the  magnitude  of  the  angle  and  how  well  the  GEOSAT  and 
buoy  wind  speed  compared.  For  angles  greater  than  1.1  degrees,  the  correlation  was 
poor.  Angles  less  than  0.4  degrees  off  nadir  exhibited  good  correlation  in  the 
GEOSAT/buoy   wind   speed  comparisons.      In  this   study,   the   nadir  angle   was  not 
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available  as  the  GEOSAT  data  used  was  preprocessed  by  Naval  Ocean  Research  and 

Development  Activity  (NORDA),  however  all  data  with  off  nadir  angles  greater  than 

f  is  not  transmitted  to  FNOC. 

B.   SSM/I  WIND  SPEED  GROUND  TRUTH  STUDY 

The  SSM/I  algorithms  have  undergone  four  revisions  since  the  system  was  put 

into  space  in  late  1987.   The  study  completed  by  Hollinger  and  Poe  [Ref.  18]  includes 

600  SSM/I  and  buoy  data  sets.    Comparisons  of  wind  speed  measurements  were  made 

within  30  minutes  and  25  kilometers  of  each  other  (Fig.  5). 
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Figure  5  SSM1  wind  speed  comparison  with  ground  truth  from  moored 
buoys  (Hollinger  and  Poe). 
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The  SSM/I  winds  are  extracted  from  a  data  set  by  the  "D-matrix"  method.  This 
is  an  empirical/statistical  approach  that  chooses  the  most  probable  atmospheric  and 
surface  conditions  that  would  produce  a  characteristic  set  of  measured  brightness 
temperatures.  The  "D-matrix"  contains  the  reference  brightness  temperatures  for  the 
SSM/I.  The  SSM/I  retrieved  brightness  temperatures  are  compared  to  the  "D-matrix" 
and  from  this  comparison  the  ocean  and  atmospheric  conditions  that  were  observed 
by  the  SSM/I  are  determined.  The  overall  scatter  of  data  points  represents  a  +\-  2 
m/s  standard  deviation  between  3  and  15  m/s.  At  this  time,  there  is  insufficient  data 
for  wind  speeds  greater  than  15  m/s  to  make  an  accurate  comparison  of  SSM/I  and 
buoy  wind  speeds. 

Interfering  factors  that  could  reduce  the  accuracy  of  the  analysis  are  heavy 
concentrations  of  water  vapor  and  the  presence  of  cloud  water  and  rain.  In  this  study, 
areas  exhibiting  these  conditions  were  included  in  the  data  to  observe  their  influence 
on  retrieved  wind  speeds. 

Positioning  error  was  discovered  by  as  much  as  20  to  30  km  for  shift  in  pixel 
latitude  and  longitude  coordinates.  However,  Hollinger  and  Poe  state  that  this  did  not 
interfere  with  the  overall  validation  of  the  wind  speed  algorithm. 
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VI.  SSM/I  AND  GEOSAT  ALTIMETER  WIND  SPEED  COMPARISON 

A.      DATA  ACQUISITION 

The  data  base  used  was  acquired  from  Fleet  Numerical  Oceanography  Center 
(FNOC)  in  Monterey,  California.  All  data  obtained  for  this  study  was  gathered  in 
near-real  time  using  a  personal  computer  (PC)  with  a  1200  baud  modem.  Co-located 
passes  of  the  DSMP  and  GEOSAT  satellites  within  one  hour  of  each  other  were 
required.  The  position  and  time  of  ascending  node  for  the  orbits  of  each  satellite 
were  calculated  well  in  advance.  These  orbit  tables  were  updated  every  two  weeks  to 
account  for  errors  in  calculating  the  exact  times  of  ascending  node  crossings. 
Adjustments  had  to  be  made  to  the  GEOSAT  orbit  table  to  account  for  the  periodic 
maneuvering  of  the  satellite  that  would  keep  it  in  a  17  day  repeat  orbit  pattern.  To 
avoid  any  microwave  interference  from  land  stations,  areas  over  the  open  sea  were 
chosen. 

Several  problems  were  encountered  in  obtaining  data  in  this  manner.  The  average 
length  of  time  to  download  a  single  SSM/I  data  set  was  at  1200  baud  was 
approximately  0.5  to  1  hour,  depending  on  its  size.  Near  real  time  SSM/I  and 
GEOSAT  data  at  FNOC  is  accessible  for  6  hours  or  less  before  it  is  removed  from 
memory.  As  a  result,  many  data  sets  obtained  were  incomplete  as  the  data  had  been 
removed  during  the  downloading  process.  I  was  unable  to  determine  at  what  times  a 
particular  orbit  would  be  resident  at  FNOC  and  for  what  length  of  time  it  would  be 
resident.    I  consequently  worked  with  the  data  available  at  the  time.    Another  difficulty 
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arose  in  the  form  of  the  GEOSAT  data  format.  The  GEOSAT  data  obtained  from 
FNOC  was  pre-processed  by  NORDA  in  Bay  St.  Louis,  Mississippi.  The  format  of 
the  GEOSAT  data  base  included  only  the  latitude,  longitude,  wind  speed,  and 
significant  wave  height.  Raw  data  would  have  been  ideal,  so  that  various  wind  speed 
algorithms  could  have  been  tested.  Occasionally  incomplete  data  sets  were  obtained 
where  wind  speeds  were  missing  for  large  areas  of  an  orbit.  This  was  due  to  data 
transmission   problems    occurring    between    NORDA    and    FNOC.  Despite    these 

difficulties,  five  data  sets  were  obtained  from  three  different  regions  of  the  globe:  the 
North  Atlantic,  the  Gulf  of  Alaska,  and  a  region  in  the  North  Central  Pacific. 

B.       DATA  FORMATS  AND  PROCESSING 

This  section  describes  the  various  formats  in  which  the  SSM/I  and  GEOSAT 
data  are  found.  In  addition,  the  process  in  which  raw  SSM/I  information  is  decoded 
and  then  combined  with  the  preprocessed  GEOSAT  wind  speeds  to  form  the  final  data 
base  will  be  discussed. 

1.      SSM/I  Data  Base 

There  are  four  forms  of  data  which  are  on  file  for  each  orbit  made  by  the 
DSMP  satellite.  All  are  encoded  in  hexadecimal  format  and  must  be  decoded  to 
decimal  before  processing. 

The  Orbit  Summary  File  (SUM)  contains  information  on  the  latitude  and 
longitude  of  each  end  of  the  sensor  swath  for  entire  satellite  orbit.  Included  with  this 
is  a  history  of  sensor  health,  i.e.,  whether  or  not  any  of  the  sensors  on  board 
experienced  any  unusual  fluxes  in  radiation,  including  the  hot  calibration  standard, 
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reference  voltages,  and  the  automatic  gain  control  From  this  file,  the  exact  time  during 
which  the  DMSP  orbit  would  coincide  with  a  GEOSAT  crossing  could  be  determined. 
This  fixed  the  amount  of  data  needed  to  be  downloaded  to  cover  as  much  of  the 
GEOSAT  track  as  possible. 

Raw  brightness  temperatures  for  each  sensor,  are  initially  stored  in  a  file  known 
as  the  Sensor  Data  Record  (SDR).  An  Antenna  Pattern  Correction  (APC)  is  applied 
to  the  brightness  temperatures  to  correct  for  spurious  energy  received  in  the  side  lobes. 
Corrected  brightness  temperatures  are  then  stored  in  the  Temperature  Data  Record 
(TDR).  It  is  in  this  format  that  brightness  temperatures  can  be  used  to  derive 
geophysical  data  such  as  precipitation  rates,  ice  edge,  ice  age,  water  vapor,  and  wind 
speed.  The  geophysical  data  are  then  stored  in  the  Environmental  Data  Record  (EDR). 
This  is  the  level  of  data  used  in  this  study. 

2.      GEOSAT  Data  Base 

GEOSAT  altimeter  data  is  stored  on  board  the  satellite  for  approximately 
12  hours  and  then  transmitted  to  a  receiving  station  at  the  Johns  Hopkins  Applied 
Physics  Laboratory  (APL)  where  it  is  archived  on  tape.  Three  data  products  are 
produced  from  the  raw  data;  The  Sensor  Data  Record,  the  Naval  Ocean  Research  and 
Development  Activity  (NORDA)  Data  Record,  and  the  Waveform  Data  Record.  The 
first  two  contain  information  about  the  measured  altitude,  significant  wave  height,  wind 
speed,  automatic  gain  control,  and  corrections  for  satellite  and  instniment  errors.  The 
last  record  contains  waveform  data  used  primarily  for  ice  studies. 

FNOC  stores  GEOSAT  data  received  from  NORDA  in  half  hour  blocks 
which  contain  latitude,  longitude,  wind   speed,   and  sea  surface  height   information. 
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The  data  is  pre-processed  and  does  not  contain  the  raw  data  necessary  to  verify  the 
algorithms  used  to  produce  the  wind  speeds  and  corrections  made  to  the  data  set. 
Reliance  was  placed  on  the  ground  truth  comparison  made  with  the  GEOSAT  wind 
speed  calculations  to  establish  the  validity  of  the  wind  speed  data. 

3.      Processing  GEOSAT  Altimeter  and  SSM/I  data  sets 

Several  Fortran  programs  were  used  to  convert  the  hexadecimal  format  of 
the  SSM/I  data  records  to  decimal.  A  program  called  COMPARE  was  written  to 
combine  the  measurements  from  both  the  GEOSAT  and  SSM/I  EDR  data  bases  for 
each  co-located  satellite  pass.  SSM/I  wind  speeds  measured  within  a  0.25  degree 
radius  of  a  GEOSAT  data  point  were  averaged  together  and  a  mean  value  and  standard 
deviation  were  obtained.  A  new  data  set  was  then  produced  which  contained  the 
GEOSAT  wind  speeds  and  the  mean  values  of  the  SSM/I  wind  analysis  from  which 
correlation  studies  were  made.  A  graphical  representation  of  the  two  data  sets  were 
produced.  The  SSM/I  EDR  also  supplied  the  rain,  cloud  water,  and  water  vapor 
analysis  necessary  to  determine  the  environmental  conditions  that  existed  during  the 
wind  speed  measurements. 

C.       DATA  ANALYSIS  AND  RESULTS 

A  comparison  of  all  the  compiled  data  sets  for  the  SSM/I  and  GEOSAT  altimeter 
wind  speed  analysis  is  seen  in  Fig.  6.  The  correlation  for  the  compiled  data  sets  is 
0.703.  From  this  graph,  four  distinct  regions  are  seen.  The  region  "A"  extends  from 
2.5  to  9  m/s.  Here,  the  SSM/I  demonstrated  good  correlation  with  the  GEOSAT 
altimeter.     Wind  speeds  between   5   and  8   m/s  represent  compiled  data  from  four 
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different  orbits,  however,  wind  speed  from  2.5  to  5  m/s  was  derived  from  a  single  data 
set,  which  resulted  in  a  region  of  sparse  data.  The  SSM/I  is  limited  to  a  minimum 
detectible  wind  speed  of  3  m/s.  This  limitation  is  due  to  the  lack  of  sensitivity  to 
changes  in  emissivity  for  wind  speeds  between  0  and  3  m/s.  The  region  greater  than 
8  m/s  demonstrated  very  poor  correlation  between  the  SSM/I  and  the  GEOSAT 
altimeter.  This  area  is  divided  in  to  three  separate  regimes.  The  region  "B" 
corresponds  to  the  bulk  of  data  which  is  greater  than  9  m/s.  Here,  the  correlation 
appears  to  be  good,  but  the  SSM/I  data  are  notably  biased  to  the  higher  wind  speeds. 
On  either  side  of  region  "B",  there  are  outlying  points.  Region  "C"  is  identified  by 
the  extremely  high  values  of  wind  speed  of  the  SSM/I  for  GEOSAT  winds  from  8  to 
12  m/s.  Region  "D"  is  represented  by  the  high  values  of  the  GEOSAT  wind  speeds 
relative  to  the  SSM/I  values  from  1 1  to  14  m/s. 

As  both  systems  have  successfully  verified  wind  speeds  in  the  regime  between 
8  to  15  m/s,  the  regions  where  the  two  measurements  differed  were  examined  in  detail. 
The  water  vapor,  cloud  water  content,  and  rain  measured  by  SSM/I  were  used  to 
determine  if  any  interference  due  to  these  substances  may  have  caused  the  bad 
correlation  in  the  high  wind  speed  areas.  This  was  done  for  each  individual  data  sets 
to  determine  if  regional  differences  may  have  made  a  contribution.  Data  sets  are 
identified  by  the  corresponding  SSM/I  orbit  number,  latitude  and  longitude  of  the 
region,  as  well  as  the  date  and  time  difference  between  the  SSM/I  and  GEOSAT  co- 
located  orbits. 
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SSM/I  VS  ALTIMETER  WIND  SPEED 
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Figure  6  SSMI  versus  GEOSAT  wind  speeds  in  meters  per  second. 
Line  number  1  represents  the  actual  correlation  of  0.703  and  line  2 
represents  a  reference  correlation  of  1.0. 

1.      Orbit  9138,  27  March,  1989 

The  scan  area  of  the  SSM/I  encompassed  a  region  bounded  by  the  following 


coordinates: 

Latitude 

Longitude 

32.50°N 

226.00°E 

32.50°N 

242.20°E 

17.00°N 

226.00°E 

17.00°N 

242.20°E 
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The  time  difference  between  the  GEOSAT  and  SSM/I  passes  is  no  greater  than  0.61 
hours.  The  track  of  the  GEOSAT  through  the  SSM/I  swath  begins  at  20°N,231.8°E 
and  ends  at  29.5°N,227.4°E.  The  wind  speeds  recorded  in  this  region  ranged  from 
2.8  to  7.9  m/s.  The  correlation  for  this  pass  is  0.925  (Fig.  7).  The  difference  in 
wind  speed  for  the  individual  satellites  demonstrates  a  maximum  of  1.9  m/s  (Fig.  8). 
The  wind  field,  as  observed  by  the  SSM/I,  is  depicted  in  Fig.  9.  The  entire  area  of 
co-located  orbits  reveals  a  region  of  high  water  vapor  content  of  approximately  60 
kg/m2  (Fig.  10).  No  cloud  water  or  precipitation  is  present  in  this  region.  Since  the 
correlation  of  the  wind  speed  for  orbit  9138  is  so  high,  one  could  safely  deduce  that 
a  water  vapor  content  as  high  as  60  kg/m2  does  not  adversely  affect  the  remote  sensing 
of  wind  speeds  by  microwave  instruments. 

2.      Orbit  8934,  13  March,  1989 

The  scan  area  of  the  SSM/I  encompassed  a  region  bounded  by  the  following 
coordinates: 

Latitude  Longitude 

57.94°N  214.33°E 

54.44°N  192.55°E 

47.39°N  216.97°E 

44.57°N  199.38°E 

The  time  difference  between  the  GEOSAT  and  SSM/I  passes  is  within  1.01  hours. 
The  track  of  co-located  orbits  begins  at  39.5°N,217.8°E  and  ends  at  27.8°N,211.5°E. 
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OOBIT    9-139 
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Figure  7  SSM/I  vs  altimeter  wind  speed.  Line  number  1  represents  the 
actual  correlation  of  0.925  and  line  2  represents  a  reference  correlation 
of  1.0. 
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Figure  8  SSM/I  and  GEOSAT  wind  speed  vs.  latitude:  orbit  91 38. 


37 


Figure  9  SSM/I  wind  speed  field  (m/s);  orbit  9138.  Line  through  field  represents 
GEOSAT  flight  path. 
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Figure  10  SSM/I  water  vapor  field  (kg/m2).  orbit  Q138.  The  line  through  the  field 
represents  the  GEOSAT  flight  path.  A  high  correlation  of  wind  speeds  exist 
between  the  GEOSAT  altimeter  and  SSM/I  despite  high  water  vapor  concentration. 
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A  break  in  the  track  was  due  to  an  interruption  in  communications  during  the  transfer 
of  data  over  a  modem.  Data  were  consequently  lost  between  31°N  and  32.5°N 
latitudes.  The  statistics  for  this  orbit  revealed  a  correlation  between  the  GEOSAT 
altimeter  and  SSM/I  wind  speeds  of  0.616.  This  data  set  spans  a  range  of  wind  speeds 
from  5  to  9.5  m/s  (Fig.  11).  The  data  points  appear  to  be  slightly  biased  toward  the 
higher  wind  speeds  for  the  GEOSAT  altimeter.  The  plot  of  wind  speed  versus  latitude 
(Fig.  12)  reveals  that  both  satellites  follow  the  trend  of  increasing  and  decreasing  wind 
speeds.  The  largest  difference  in  wind  speeds  (approximately  2.5  m/s)  appears  at 
29.2°N,  32.5°N,  35.5°N,  and  38°N  latitudes.  The  images  of  wind  speed  (Fig.  13), 
water  vapor  (Fig.  14),  and  precipitation  (Fig.  15)  show  that  the  track  of  the  co-located 
orbit  does  not  pass  through  any  interfering  atmospheric  substances  and  that  there  were 
no  unusual  gradients  in  the  vicinity.  Overall,  orbit  8934  demonstrated  good  correlation 
between  the  GEOSAT  altimeter  and  the  SSM/I. 


3.      Orbit  8058(7). 

,  10  January,  1989 

The  scan  area  < 

:>f  the  SSM/I  encompassed  a  region  bounded  by  the  following 

coordinates: 

Latitude 

Longitude 

59.73°N 

217.54°E 

57.84°N 

197.63°E 

51.88°N 

223.64°E 

48.89°N 

205.01  °E 
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Figure  11  SSM/I  vs.  GEOSAT  WIND  SPEED:  ORBIT 
8934.  Line  1  represents  the  linear  regression  line.  Line 
2  is  a  reference  slope  of  1.0. 
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Figure  12  SSM/I  and  GEOSAT  winds  vs.  latitude: 
8934. 


orbit 
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Figure  13     SSM/I  wind  speed  field  (m/s):  orbit  8934. 
field  represents  the  GEOSAT  flight  path. 


Line  through  the  wind 
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Figure  14      SSM/I   water  vapor  field   (kg/m2);  orbit  8934. 
represents  the  GEOSAT  flight  path. 


Line  through   field 
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Ficure  15     SSM/I  precipitation  field  (mm/hr);  orhit  8934.     Line  through  field 
represents  the  GEOSAT  flight  path. 
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The  time  difference  between  the  GEOSAT  and  SSM/T  passes  is  no  greater  than  1.10 
hours.  The  track  of  the  GEOSAT  through  the  SSM/I  swath  begins  at  49°N,219°E  and 
altimeter  for  this  particular  orbit.  By  inspection,  the  correlation  between  the  two 
instruments  is  poor.  The  calculated  correlation  is  0.208.  Fig.  17  illustrates  the 
variation  in  wind  speeds  going  from  latitude  49°N  to  57.5°N.  Beginning  at  the  most 
southern  point,  the  difference  between  the  SSM/I  and  GEOSAT  altimeter  grow  from 
approximately  1  m/s  to  as  much  as  3.5  m/s.  The  corresponding  SSM/I  wind  field  is 
represented  in  fig.  18.  Examination  of  the  SSM/I  water  vapor  field  (Fig.  19),  shows 
the  orbit  tracks  cross  through  a  gradient  of  water  vapor  which  changes  from 
approximately  10  to  30  kg/m2  over  a  distance  of  about  50  nm.  The  track  continues 
through  an  area  of  saturated  water  vapor,  i.e.,  greater  than  60  kg/m2  and  then  exits  into 
a  region  of  approximately  30  kg/m2.  The  crossing  of  the  gradient  at  51°N  coincides 
with  the  rapid  increase  in  the  SSM/I  wind  speed,  while  the  altimeter  remains  relatively 
constant.  Examination  of  the  SSM/I  precipitation  field  (Fig.  20)  reveals  precipitation 
rates  from  1  to  4  mm/hr  at  51°N  and  from  53.5°N  to  56°N,  along  the  co-located  orbit 
track.  From  56°N  to  the  end  of  the  track,  cloud  water  is  present  up  to  0.125  kg/m2 
(Fig.  21).  This  suggests  that  in  regions  of  cloud  water  and  precipitation,  SSM/I  winds 
are  greater  than  GEOSAT  values. 
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Figure  16  SSM/I  vs.  GEOSAT  wind  speed:  orbit  80587.  Line  1 
represents  an  actual  correlation  of  0.208,  while  line  2  represents  the 
reference  correlation  of  1.0. 
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Figure  17  SSM/I  and  GEOSAT  wind  speed  vs.  latitude:  orbit  80587. 
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Figure  18  SSM/I  wind  field  (m/s);  orbit  8058(7).  Line  through  field  represents  the 
GEOSAT  flight  path. 
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Figure  19  SSM/T  water  vapor  field  (kg/m2):  orbit  8058(7).     Line  through  field 
represents  the  GEOSAT  flight  path. 
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Figure  20  SSM/T  precipitation  field  (mm/hr);  orhit  80^8(7) 
field  represents  the  GEUSA'I   flight  path. 


Line  through  wind 
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Figure  21   SSM/I  cloud  water  field  (kg/m2);  orbit  R05R(7). 
represents  the  GEOSAT  flight  path. 


I  ine  through  field 
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4.      Orbit  8058(6),  10  January,  1989 

The  scan  area  of  the  SSM/I  encompassed  a  region  bounded  by  the  following 
coordinates: 

Latitude  Longitude 

51.88°N  223.64°E 

48.89°N  205.01  °E 

42.14°N  225.01°E 

39.60°N  210.1 1°E 

The  time  difference  between  the  GEOSAT  and  SSM/I  passes  is  within  1.00  hours. 
The  track  of  the  GEOSAT  through  the  SSM/I  swath  begins  at  40.8°N,224.6°E  and 
ends  at  49.5°N,218.6°E.  The  corresponding  graph  of  SSM/I  vs  GEOSAT  wind  fields 
is  presented  in  Fig.  22.  This  data  set  spans  a  range  of  wind  speeds  from  7  to  12  m/s. 
The  correlation  in  this  graph  is  good  below  8  m/s  and  then  the  SSM/I  appears  to  read 
slightly  higher  then  the  GEOSAT  altimeter.  Overall,  the  correlation  is  calculated  to 
be  0.935.  A  graph  of  wind  speed  versus  latitude  (Fig.  23)  shows  that  both  instruments 
were  within  1  m/s  from  41°N  to  about  46.5°N  and  afterwards  the  difference  began  to 
open  to  about  2  m/s  near  48.5°N.  Examination  of  the  SSM/I  wind  field  (fig.  24) 
shows  that  the  location  where  the  difference  in  wind  speed  begins  to  grow  occurs 
around  the  same  area  that  the  wind  speeds  begin  to  increase  rapidly.  The  SSM/I  water 
vapor  field  (Fig.  25)  does  not  demonstrate  any  significant  gradient  which  would 
interfere  with  the  measurement  of  wind  speed  or  any  area  on  increased  water  vapor. 
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Figure  22  SSM/I  vs.  GEOSAT  wind  speed:  orbit  80586.  Line  l 
represents  an  actual  correlation  of  0.935,  while  line  2  represents  a 
reference  correlation  of  1.0. 
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Figure  23  SSM/I  vs.  GEOSAT  winds  vs.  latitude:  orbit  8058(6). 
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Fiunre  24    SSM/T   wind   speed   field   (m/s):    orbit    R058(fi). 
GEUSAT  flight  path. 


line   represents  the 
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Fimire  25   SSM/I  water  vapor  fieM  (kg/m2):  orbit  8058(6). 
represents  the  GEOSAT  flight  path. 


I  ine  through  field 
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Figure  26  SSM/I  precipitation  fiHrl  (mtn/hr):  orhit  80^8(6).  The  line  through  the 
field  represents  the  GEOSAT  flight  path. 
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There  is,  however,  a  gradual  increase  in  water  vapor  content  in  the  northeast  quadrant 
of  the  field  from  15  to  30  kg/m2.  The  precipitation  field  (Fig.  26)  shows  precipitation 
to  the  northeast  of  the  GEOSAT-SSM/I  track,  but  again,  there  is  no  evidence  of 
precipitation  occurring  along  the  track.  Overall,  there  appeared  to  be  no  interfering 
atmospheric  substances  in  the  region  of  the  co-located  orbits. 

5.      Orbit  7925,  2  January,  1989 

The  scan  area  of  the  SSM/I  encompassed  a  region  bounded  by  the  following 
coordinates: 

Latitude  Longitude 

79.98°N  335.99°E 

70.99°N  296.68°E 

52.27°N  348.24°E 

49.25°N  329.42°E 

The  time  difference  between  the  GEOSAT  and  SSM/I  passes  is  less  than  0.50  hours. 
The  track  of  the  GEOSAT  through  the  SSM/I  swath  begins  at  54.6°N,343°E  and 
terminates  at  65.4°N,324.4°E.  The  plot  of  SSM/I  versus  GEOSAT  wind  speeds  (Fig. 
27)  demonstrates  a  tremendous  amount  of  scatter  in  the  data.  Wind  speeds  range  from 
7.5  m/s  to  20.5  m/s.  The  correlation  calculated  for  this  data  set  is  -0.301.  Looking  at 
the  plot  of  wind  speed  versus  latitude  (Fig.  28),  one  can  see  that  there  are  two  distinct 
regions  where  the  difference  in  wind  speeds  are  greater  than  3  m/s.  This  occurs 
between  57.5°N  and  60°N  and  again  from  63.5°N  until  65.3°N.  The  SSM/I  wind  field 
(Fig.  29)  exhibits  a  range  of  wind  speeds  from  5  to  31  m/s.    The  highest  wind  speeds 


56 


occur  along  a  band  extending  from  the  south  central  region  of  the  field  toward  the 
coast  of  iceland.  High  wind  speeds  are  also  observed  along  the  coast  of  greenland. 
The  SSM/I  water  vapor  field  (Fig.  30)  exhibits  a  high  concentration  of  water  vapor  in 
the  southeast  quadrant  of  the  SSM/I  and  GEOSAT  orbits.  The  large  white  area  in  the 
water  vapor  image  represents  values  higher  than  60  kg/m2.  An  examination  of  the 
precipitation  field  (Fig.  31)  discloses  a  major  region  of  precipitation  occurring  along 
the  track  as  well.  The  location  of  the  precipitation  coincides  with  the  same  area  of 
Fig.  28  that  reveals  a  large  difference  in  wind  speed  between  the  SSM/I  and  GEOSAT 
altimeter  (57.5°N  to  60.0°N).  The  extent  of  precipitation  appears  to  cover  too  small 
a  region  to  be  considered  the  only  reason  for  the  disparity  between  the  GEOSAT  and 
SSM/I  wind  speeds.  The  SSM/I  cloud  water  field  (Fig.  32)  shows  that  this  entire 
region  is  covered  by  high  cloud  water.  From  this  evidence,  the  cloud  water  and 
precipitation  together  have  apparently  caused  the  large  difference  in  wind  speed  values 
between  the  SSM/I  and  GEOSAT  altimeter. 

A  second  area  of  concern  lies  between  63.5°N  and  65.3°N  (Fig.  28).  This  area 
is  very  close  to  the  coast  of  Greenland,  approximately  120  km  at  its  nearest  point. 
The  altimeter  wind  speed  measurement  steadily  increases  as  it  approaches  Greenland, 
while  the  SSM/I  wind  speed  measurements  decrease.  This  phenomena  may  be  due  to 
the  closeness  of  the  land  mass  and  side  lobe  effects  on  the  SSM/I.  Criteria  for 
nearness  to  land  is  not  clear  at  this  time.  Chelton  and  McCabe  [Ref.  18]  and  Wentz 
et  al  [Ref.  19]  consider  only  those  regions  where  microwave  remote  sensing  is  done 
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Figure  27  SSM/I  vs.  GEOSAT  wind  speed:  orbit  7925.  Line  1 
represents  an  actual  correlation  of  -0.301,  while  line  2  represents  a 
reference  correlation  of  1.0. 
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Figure  28  SSM/I  and  GEOSAT  wind  speeds  vs.  latitude. 
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Figure  29      SSM/T   wind   zpeed   field   (m/s):    orbit    7925. 
GEOSAT  flight  path. 


Line   represents   the 
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Figure  30   SSM/I   water   vapor   field    (kg/m2);    orbit   7025.      I,ine   through    field 
represents  the  GEOSAT  flight  path. 
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Fimire  31    SSM/1   precipitation   field   (mm/hr);   orbit   7925.      Line  through   field 
represents  GEOSAT  flight  path. 
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Figure  32    SSM/I    cloud    water   field    (kg/m2):    orhit   7Q25. 
represents  the  GEOSAT  flight  path. 


Line    through    field 
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from  200  to  800  km  from  the  nearest  land  mass.  They  sight  erratic  behavior  due  to 
side  lobes  and  the  change  in  emissivity,  when  making  microwave  measurements  of  the 
open  ocean  are  made  close  to  land. 
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VII.  CONCLUSIONS 


To  assess  the  relationship  between  the  wind  speeds  as  measured  by  the 
GEOSAT  altimeter  and  the  SSM/I,  the  correlation  (R),  R2,  slope  and  intercept  were 
calculated.    Table  one  contains  the  results  of  the  analysis  of  each  data  set. 

The  value  of  R2  represents  the  proportion  of  total  variance  in  the  values  of  the 
SSM/I  wind  speed  data  that  can  be  accounted  for  by  a  linear  relationship  to  the 
GEOSAT  wind  speeds.  In  orbits  9138  and  8058(6),  more  than  90%  of  the  variance 
is  explained  by  this  linear  relationship.  Orbits  7925  and  8058(7)  demonstrate  almost 
no  correlation  between  the  SSM/I  and  altimeter  wind  speeds. 

The  intercept  gives  the  bias  between  wind  speed  measurements.     Orbit  9138 

and  8058(6)  exhibit  the  smallest  bias,  approximately  1  m/s.    The  accuracy  goal  of 

Table  1.  Each  data  set  is  identified  by  its  orbit  number  along  with  its 
associated  correlation  (R),  R2,  intercept,  slope,  and  atmospheric  parameters 
present  during  the  wind  speed  measurement. 


ORBIT 

R 

E2 

SLOPE 

INTERCEPT 

COMMENTS 

9138 

0.925 

0.8583 

1.284 

-1.297 

High  water  vapor 
concentration 

8934 

0.616 

0.5070 

0.644 

6.889 

Clear  Atmosphere 

8058(7) 

0.208 

0.0431 

0.329 

10.599 

Precipitation,  cloud  water,  and 
water  vapor 

8058(6) 

0.935 

0.8739 

1.247 

-1.800 

Clear  Atmosphere 

7925 

0.301 

0.0900 

0.258 

16.980 

Precipitation,  cloud  water,  and 
water  vapor 
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the  GEOSAT  altimeter  and  the  SSM/I  has  been  2  m/s.  Both  orbits  9138  and 
8058(6)  are  within  that  goal.  Orbits  7925,  8058(7),  and  8934  have  a  bias  in  wind 
speed  greater  than  2  m/s.  The  slope  is  also  a  good  indicator  of  the  relationship 
between  the  SSM/I  and  GEOSAT  altimeter  wind  speeds.  A  perfect  relationship 
between  wind  speed  measurements  would  give  a  slope  of  1.0.  Orbits  9138  and 
8058(6)  have  a  slope  very  close  to  one,  while  7925  and  8058(7)  have  a  slope  near 
0.0. 

There  are  four  conclusions  that  are  made  from  the  analysed  data  in  this  paper: 

(1)  The  data  sets  which  did  not  contain  cloud  water  or  precipitation  along  the 
co-located  SSM/I  and  GEOSAT  orbits  (orbits  9138,  8934,  and  8058(6)), 
demonstrated  a  good  correlation  between  the  two  systems. 

(2)  Although  a  large  amount  of  water  vapor  was  present  along  the  entire  track 
of  orbit  9138,  there  appeared  to  be  no  evidence  of  bias  with  the  measured  wind 
speeds  from  either  instrument.  Since  the  SSM/I  corrects  for  the  presence  of  water 
vapor  in  the  atmosphere,  it  is  obvious,  from  the  high  correlation  between  the 
GEOSAT  altimeter  and  the  SSM/I,  that  the  SSM/I  water  vapor  correction  works 
well.  The  altimeter  was  not  corrected  for  water  vapor,  suggesting  that  water  vapor 
corrections  on  the  altimeter  were  small. 

(3)  The  data  sets  which  contained  cloud  water  and  precipitation  (orbits  8058(7) 
and  7925),  displayed  large  differences  in  recorded  wind  speeds  between  the  SSM/I 
and  GEOSAT  altimeter.  In  these  areas,  the  SSM/I  wind  speed  would  increase, 
while  the  altimeter  wind  speed  decreased  relative  to  one  another  and.  consequently, 
very  poor  correlation  resulted.     There  are  several  explanations  for  the  difference  in 
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wind  speed  relative  to  each  other  experienced  by  both  sensors.  The  increased 
emissivity  associated  with  precipitation,  or  cloud  water,  is  the  most  likely  cause  of 
increased  wind  speeds  experienced  by  the  SSM/I,  by  introducing  a  "warm 
atmospheric  term.  The  GEOSAT  may  experience  a  decrease  in  wind  speed  at  the 
same  location.  This  may  have  been  a  result  of  a  decrease  in  ocean  surface 
roughness.  The  reduction  in  surface  roughness  is  caused  by  rain  drops  impacting  on 
the  ocean  surface  and  consequently  attenuating  the  surface  gravity  waves  [Ref.  21]  . 
The  reduction  in  surface  roughness  would  result  in  a  larger  backscatter  cross  section 
for  the  GEOSAT  altimeter.  This  reduced  backscatter  is  interpreted  as  a  decrease  in 
wind  speed,  which  is  what  is  observed  in  the  area  of  precipitation  in  orbits  7925 
and  8058(7).  An  alternative  explanation  for  the  low  GEOSAT  wind  speed  in  the 
area  of  precipitation  is  that  the  measurement  was  correct  and  that  the  SSM/I  had  a 
bias  to  high  wind  speeds  in  the  presence  of  rain. 

(4)  Erratic  behavior  is  seen  in  microwave  measurements  of  the  open  ocean 
near  a  land  mass.  This  is  evident  in  orbit  7925  when  measurements  of  wind  speed 
and  atmospheric  properties  were  made  within  120  km  of  the  Greenland  coastline. 
The  resulting  data  from  the  coastal  areas  was  highly  erratic.  Based  on  only  one 
data  set  near  a  coastal  feature,  it  appears  that  wind  speed  measured  within  200  km 
of  a  land  mass  is  not  reliable.  This  is  most  likely  due  to  the  radiance  detected  by 
antenna  sidelobes  which  mask  the  actual  measurement  of  radiation  from  the  main 
lobe  of  the  antenna. 

The  wind  speed  verification  studies  of  the  GEOSAT  altimeter  and  the  SSM/I 
demonstrate  that  both  active  and  passive  microwave  systems  measure  similar  values 
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of  wind  speed  over  the  ocean.  These  instruments  are  comparable  in  open  ocean 
regions  at  least  200  km  from  land  and  free  of  cloud  water  and  precipitation. 
Utilization  of  GEOSAT  altimeter  and  SSM/1  wind  speed  data  in  a  numerical 
forecasting  model  requires  that  a  system  for  screening  wind  data  for  cloud  water, 
precipitation,  and  proximity  to  land  masses  be  used. 
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